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SMALL  SIGNAL  ANALYSIS  OF  THE  INDUCED 
RESONANCE  ELECTRON  CYCLOTRON  MASER 


I.  Introduction 

Generation  of  intense  radiation  in  the  aicrovave  regime  utilizing 

electron  cyclotron  interaction  has  been  proposed  independently  by  a 

1-4 

number  of  researchers  in  the  late  1950's.  Electrons  gyrating  in 
resonance  vith  the  radiation  field  can  experience  a  bunching  in  the 
relative  vave-particle  phase  through  the  dependence  of  the  cyclotron 
frequency  on  the  relativistic  mass.  High  amplification  of  the 
radiation  field,  known  as  easing  action,  results  for  Doppler  shifted 
frequencies  slightly  above  the  electron  cyclotron  frequency.  Electron 
cyclotron  masers,  also  called  gyrotrons,^"^  have  demonstrated 
efficient  high  pover  generation  of  electromagnetic  vaves  at  centimeter 
wavelengths. 

For  many  purposes  it  is  of  practical  interest  to  develop  high 
power  generation  capability  at  millimeter  and  submillimeter 
wavelengths.  Potential  areas  of  application  include  advanced 
accelerators,  short  wavelength  radar,  plasma  heating  in  fusion 
reactors  and  spectroscopy.  The  shortest  vavelength  for  single  mode 
operation  in  a  closed  resonator  is  tied  to  the  transverse  dimension  of 
the  cavity.  For  radiation  vavelengths  much  shorter  than  the 
transverse  dimensions,  a  multimode  excitation  will  result  from  the 
small  frequency  separation  among  nearby  modes.  The  node  selectivity 

is  greatly  improved  by  the  use  of  an  open  resonator  configuration,  the 

.  i  19, 20 

quasi-optical  maser. 

A  new  configuration  has  recently  been  proposed‘^’^  which 
utilizes  the  benefits  of  the  open  resonators  and  at  the  same  time 
minimizes  the  detrimental  effects  of  the  injected  electron  beam  energy 
spread.  The  operating  frequency  in  the  induced  resonance  electron 
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cyclotron  (IREC)  quasi-optical  maser  is  upshifted  by  a  factor  y 

relative  to  the  relativistic  electron  cyclotron  frequency.  It  has 

been  shown  that  for  operation  at  the  optimum  index  of  refraction  the 

efficiency  is  relatively  insensitive  to  the  beam  energy  spread  and 

the  sensitivity  to  the  effect  of  pitch  angle  spread  can  be  minimized. 

The  index  of  refraction  is  adjustable  by  varying  the  angle  between  the 

resonators  (see  Fig.  (1))  and  the  guide  field,  and  can  be  chosen  to 

minimize  the  effects  of  finite  beam  quality.  Finally,  by  spatially 

tapering  the  magnetic  field  the  operating  efficiency  can  be  increased. 

In  this  paper  we  limit  ourselves  to  analyzing  the  small  signal 

efficiency  characteristics  of  such  a  device.  Ve  include  the  effects 

of  the  Gaussian  profile  for  the  radiation  envelope  considering  a 

uniform  magnetic  field  for  simplicity.  Nonlinear  effects  and  the  role 

30 

of  the  magnetic  field  tapering  are  treated  elsewhere. 

The  remainder  of  this  paper  is  organized  as  follows.  In  Sec.  II 
we  describe  the  field  configuration  and  the  equations  of  motion.  In 
Sec.  Ill  we  derive  the  linear  energy,  power  efficiency  and  start-up 
current  condition.  In  Sec.  IV  numerical  results  and  conclusions  are 
presented . 


II .  Field  and  Particle  Dynamics 

The  configuration  for  the  induced  resonance  electron  cyclotron 
(IREC)  quasi-optical  maser  is  shown  schematically  in  Fig.  1.  The 
interaction  cavity  is  formed  by  two  quasi-optical  resonators 
intersecting  at  an  angle  2a  where  a  is  the  angle  relative  to  the 
external  magnetic  field  Bq  in  the  z  direction. 


The  bean  radius  is  much  smaller  than  the  Gaussian  width  rQ  (spot 
size)  for  the  radiation  envelope.  In  the  limit  of  small  Larmor  radius 
p  compared  to  the  perpendicular  vavelength  k  3  «  1  ve  can  approximate 
the  vector  potential  in  the  interaction  regime  by 

A.J.  -  A^z)  exp[i*(z,  t)]  \  (ex  +  iey)  +  cc 

♦  al(z>  exp[ i*(z, t) ]  j  (ex  -  iey)  +  cc.  (1) 

Since  ve  are  interested  in  the  synchronous  interaction  of  the  gyrating 
electrons  with  the  radiation,  ve  have  kept  only  the  forward 
propagating  wave  component  *(z,t)  «  kzz  -  <ot  ♦  ♦  .  The  amplitudes  A^ 
and  A^  for  the  right-  and  left-handed  polarized  wave  component, 
respectively  are  given  by 

Aj^, L^2)  m  A^expl-  z2/L2), 

^r.L  “  Ao  <COSa  *  1)’ 

L  »  r  /sina,  (2) 

o 

where  Aq  and  rQ  are  the  amplitude  and  spot  size  for  each  individual 
resonator  beam. 

Ve  use  the  guiding  center  description  for  the  particle  orbits 

x  »  x  ♦  p  sinC,  y  *  y  -  p  cost, 

g  S 

p  «  p  ♦  p.cosC,  p  *  p  ♦  p.sinC,  (3) 

Kgx  yL  Ky  *gy 

to  obtain  the  nonlinear  relativistic  equations  of  motion.  In  this 

representation  (x  ,y  )  and  (p  ,p  )  denote  the  transverse  coordinates 

K  g  g  g*  Hgy 


and  momentum  of  the  particle's  guiding  center,  p  is  the  Larmor  radius, 

p^  is  the  aagnitude  of  the  transverse  momentum  and  £  is  the  momentum 

space  angle.  Ve  assume  that  x,  y,  p  ,  p  ,  p  and  p,  are  slowly 

x  y  x 

changing,  on  the  spatial  scale  of  a  gyroperiod.  An  additional 
condition  for  ignoring  finite  k^  effects  is  that  the  guiding  center 
shift  in  the  x  direction  be  small  k^Ax  <<  1,  valid  for  a  «  1  where  k^ 
is  k  sina.  Using  the  Lorentz  force  equation  together  with  Maxwell's 
equations  and  retaining  only  the  right-hand  polarized  wave  component 
the  nonlinear  relativistic  equations  of  motion  are  cast  into  the  form 

u.  '  -  -  l(wy/cu ,)-k _]a(z)cos*  a'(z)sin^,  (4a) 

X  z  z 

u  '  -  -  (u./u,)|k  a(z)cos*  +  a'(z)sin*],  (4b) 

Z  X  z  z 

♦'  -  -  (yAw/cuz)  ♦  (1/u^) ^(»v/cuz)-kzja(z)sin*  +  a'(z)cos*j,  (4c) 

The  prime  (')  in  Eqs.  (4)  signifies  the  d/dz  derivative,  u  « 

2  2  1/2 

p/m  c  -  yv/c,  y  *  (1  ♦  u.  ♦  u_ )  is  the  relativistic  mass  factor, 
o  •*  X  z 

2 

.<*)  -  MV  z)/mQc  is  the  normalized  radiation  amplitude,  ♦  -  C  +  ♦  is 

the  relative  phase  between  the  radiation  field  and  particle, 

n  -  ck  /«  *  cosa  is  the  refractive  index  associated  with  the  radiation 
z 

field,  Aw  -  (w(l-n0z)  -  2q/ y1  is  the  frequency  mismatch  term  and 

Qq  m  |e|Bo/moc  is  the  nonrelativist ic  electron  cyclotron  frequency.  Using 

Eqs.  (4)  the  rate  of  change  of  y  is  given  by 

y'  »  -  uKu./cu  )a(z)cos*.  (5) 

X  z 

The  frequency  mismatch  Aw  and  its  dependence  on  the  particle 
energy  through  the  relativistic  correction  y,  provide  the  mechanism 
for  the  masing  action  (phase  bunching). 


III.  Efficiency 

One  of  the  central  issues  concerning  maser  operation  is  the 

efficiency  of  the  configuration.  Efficiency  calculations  have  been 

carried  out  for  various  configurations  in  the  general  categories  of 

5-18  21-28 

the  closed  resonator  gyrotron  *  or  the  open  resonator  quasi- 

19  20 

optical  maser.  ’  Vhile  it  is  generally  recognized  that  nonlinear 

saturation  mechanisms  are  very  important  for  the  full  power  operation, 

it  is  useful  to  carry  out  the  small  signal  efficiency  calculation  in 

order  to  compute  the  start-up  current.  Expressions  for  the  small 

amplitude  efficiency,  obtained  in  closed  form,  provide  some  guidelines 

in  selecting  the  optimum  operating  parameters. 

Assuming  steady  state  operation,  with  the  number  of  particles 

crossing  the  unit  area  per  unit  time  n  v  being  constant,  the 

o  z 

efficiency  can  be  defined  by 

*e  ■  -  '  <ro  - 1>_1  Fu„<i!o>fd*  If  •  <6> 

_  00 

In  Eq.  (6),  the  bracket  <  >  signifies  the  average  over  the  initial 

distribution  in  phase  space,  the  subscript  ±®  stands  for  the  initial 

and  final  values  at  z  =  ±®  respectively  and  3y/3z  is  a  function  of  the 

initial  conditions  r  =  y(z>P^0*PZ0»^0)  •  t*le  cold  beam  limit  with 

the  initial  distribution  function  given  by  fQ< p^ , pz , =  (nQ/2»ip^) 

&(p.  -  p.  )5  (p  -  p  )  the  average  reduces  to  an  average  over 
X  10  z  zo 


a.  Small  Signal  Efficient 


Ve  proceed  to  compute  the  small  signal  power  efficiency  by 
evaluating  the  right-hand  side  of  Bq.  (6)  using  Eq.  (5).  A  first 
order  expansion  for  the  quantities  u^  =  uj/°^  +  r  =  + 

will  suffice  for  a  quadratic  expression  in  the 
wave  amplitude  a.  The  integrand  on  the  right-hand  side  of  Eq.  (6)  is 
expanded  using  the  linear  solutions  from  Eqs.  (4a)-(4c).  The 
evaluation  of  the  final  result  is  considerably  simplified  by 
performing  the  phase  space  average  over  the  angle  <^0  before  the 
spatial  integration  over  z.  Expanding  the  products  of  the 
trigonometric  terms  inside  the  integral  in  Eq.  (6)  into  sums  and 
averaging  over  leads  to 


(j dz  li)  ■  -  fey  [(» *  \  y)i  d2a<z>i d2'  sin  do<2-2,) 


ZO  -® 


(U  2  uyy  \  ”  00 

kz(l  +  j  jj  dza(z) J  dz'  a(z')cosA  (z-z') 

UZO  20  -«  -® 


i  r  2  u.  u, 

1  l«_  lo  .  lo 

"212  2  "  z  2 

vc  u„  u 

ZO 


2  ' 
O 


z  z' 


VCU  CU  J . 

ZO  ZO  ZO  - 


\  €*  C* 

ij  dza(z)J  dz'J  dz”a(z")sinAo(z-z") 


-00  _00 


Z  Z' 


+  \  p2-  -  — —If  dza(z)  T  dz'  [  dz"  da-g>  cosA  (z-z")l ,  (7) 

u  cu  cu_  /J_  J  J 

ZO  ZO  — 
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where  A 
o 


^-2.  a(z)  =  aoexp[-z2/L2] 
ZO 


Ve  evaluate  the  remaining  integrals  in  Eq.  (14)  and  express 

the  final  result  in  terms  of  the  parameters  £  -  wx  *  (wy  /cu  )L,  x 

o  zo 

being  the  transit  time  through  the  interaction  regime,  and  the 
relative  frequency  mismatch  Awq/w.  We  find 


2  -2 


2  w1’ 


& 


12@J<1  -  n2) 


(‘  *  *9]  ^  - 


e2  ng  j;2 

o  zo 


e*r 


9 


(8) 


with  0  ■  u,  /u  ,  the  initial  pitch  angle, 
o  J.O  zo 

2  2  2 

The  efficiency  is  proportional  to  exp[-  1/2  ^  Awq/w  ]  where 
exponent  £(Awo/w)  is  equal  to  Awqx,  the  advance  in  the  relative  phase 
A<J<o  between  the  wave  and  the  particle  over  the  interaction  regime. 

For  typical  values  of  £  »  1  and  Awq/w  <<  1  the  expression  in  braces 
in  Eq.  (8)  is  simplified  to 


1  + 


s2e 2  aw 

o  lo  _ o 

2  w 


(9) 


2  2  2 

where  =  £  (1-n  )  is  independent  of  a.  In  (9)  we  have  omitted  the 

small  terms  that  originate  from  the  gradient  terms  3a/3z  in  the 
equations  of  motion.  Treating  (9)  as  a  quadratic  form  in  Aw/w  we  find 
the  regime  for  positive  efficiency,  given  by 


.-1  Aw 


-1 


2(‘  -  i1  *  e0  Ko)  KoEo)  <  IT  <  9&o)  •  <10> 
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The  upper  limit  in  A«o/g)  is  due  to  a  finite  n  and  results 

2 

from  the  negative  contribution  of  the  quadratic  term  (Au/m)  that 
overtakes  the  positive  contribution  of  the  linear  term  Aco/w  for  small 
angles  sin2a  <  (2n926  /0.2)  (A<o  /«). 

In  order  to  determine  the  maximum  efficiency  within  the 
positive  regime,  we  parameterize  Eq.(8)  as  a  function  of  x  =  £  Aco/w, 
since  the  exponential  is  the  main  factor  limiting  efficiency.  Setting 
dri/dx  -  0,  we  obtain 

3  2 

CjX  -  Cj  X  -  CjX  +  cQ  =  0  ,  (11) 


with  Cj  = 

c2  *  co  “ 
Cj  and  c^ 


2  2 
(1  +  30  )0  cosa  -  1,  c,  =  0  0  cosa  and 

O  ZO  J  0  ZO 

2 

(1/2)  3jl0^0  sina.  Observing  that  the  terms  proportional  to 
can  be  omitted  provided  that  cq  =  »  c^  -  c^  -  1  or 


sina  » 


e23 

o  zo 
2  ’ 
*Io<o 


(12) 


we  can  show  that  x  =  1.  In  the  special  case  x  =  1,  we  obtain  the 
maximum  efficiency 


^ax 


it  2  -1/2 

T  a«e 
A  0 


a2  i3 
PloS 

Yo(V1) 


(sina) 


-1 


(13) 


The  overall  efficiency  increases  with  decreasing  a  (increasing  index 
of  refraction)  provided  that  inequality  Eq.  (12)  remains  valid.  For 
very  small  a  Eq.  (13)  fails  and  a  solution  of  the  cubic  Eq.  (11)  is 


necessary. 


b.  Start-up  Current 


We  are  in  position  now  to  calculate  the  start-up  beam  current 
utilizing  the  power  efficiency  coefficient.  Amplification  of  the 
electromagnetic  field  energy  will  result  if 

nPb  >  |f,  <H) 


where  £  is  the  total  electromagnetic  energy  stored  in  both  cavities 
e  =  JUDdV  =  2V(a>2/c2)(A2/4n) ,  V=nr2  L_,  de/dt  =  (w/Q )e,  Q  is  the 

K  0  0  i 

quality  factor  for  the  cavity  and  is  the  electron  beam  power. 

The  optimum  power  efficiency  hmax  is  given  by  Eq.  (13).  The 
cavity  Q  is  given  by 


Q  = 


2  n 


4 

A  ’ 


(15) 


where  1^,  is  the  effective  resonator  length  and  X  the  wavelength. 
Combining  Eqs.  (13),  (14),  (15)  and  expressing  Aq  in  terms  of  aQ  from 
Eq.  (2)  we  obtain 


P  v  X  ..  _  .  exp^  moc  0zoyo(yo  ^  2sina  .... 

Pb  >  r  (1‘Ref)  ~2—  -~2  - 3 -  - - 72’  (16) 

o  4  Jt  |  e  I  ei0  (1+cosa) 


where  ?b  = 
beam.  For 


0.1, 

angle 


PbVb’  *s  t*ie  current  anc*  Is  v°ltage  of  the  electron 
typical  parameters  =  0.25  x  lO^eV,  X/rQ  =  10-*,  1-R^  = 
1.5,  0zo  =  0.64,  0^o  =  (/3yo)-^,  and  the  optimum  operation 

45°,  the  start-up  current  is 


Ib  >  4.6  A. 
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IV.  Conclusion  and  Summar; 


Ve  have  performed  the  small  signal  analysis  for  an  oscillator 
configuration  capable  of  generating  radiation  in  the  millimeter  and 
the  submillimeter  regime.  The  threshold  for  the  start-up  current  was 
found  to  be  veil  within  the  existing  capabilities  of  today's  long 
pulse  mildly  relativistic  beams.  Our  theoretical  linear  efficiency 
results  are  plotted  as  solid  lines  in  Figs.  2-4  against  the  numerical 
results  (dots)  obtained  by  direct  integration  of  the  fully  nonlinear 
Eqs.  (4)  for  small  wave  amplitude.  Plots  of  the  linear  efficiency  as 
a  function  of  the  controlling  parameter  £  &«/«  for  constant  radiation 
amplitude  aQ  and  constant  spot  size  rQ  are  shown  in  Fig.  2,  with  each 
curve  corresponding  to  a  different  index  of  refraction  n  =  cosa.  The 
maximum  efficiency  for  all  plots  occurs  at  ^  A»/a>  =  1  in  agreement 
with  Eq.  (13).  Small  signal  efficiency  increases  with  increasing 
n  =  cosa  roughly  proportionally  to  the  length  of  the  interaction 

29-30 

regime  L  =  rQ/sina.  In  Fig.  3,  the  optimum  index  of  refraction 
2 

n  =  3  /(1-3,  ),  to  minimize  the  effects  of  beam  energy  spread,  is 

20  -LO 

held  constant,  and  the  interaction  length  L  is  changed  by  increasing 

the  width  of  the  radiation  envelope  rQ.  Figure  4  is  a  comparison  of 

the  theoretical  small  signal  efficiency  with  the  numerically 

calculated  nonlinear  efficiency  as  a  function  of  wave  amplitude  aQ. 

-4 

The  agreement  is  good  for  aQ  <  3  x  10  .  Nonlinear  saturation  occurs 

_3 

for  aQ  >  1  x  10  .  Obtaining  the  scaling  of  the  efficiency  in  the 

nonlinear  regime  is  not  possible  analytically.  Numerical  studies  of 
the  high  power  performance,  however,  have  demonstrated  good  nonlinear 
efficiency. 
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Figure  3.  Theoretical  (solid  line)  and  nuaencal  (dots)  plots  of 
linear  efficiency  n  versus  1  »  for  various  ('.aussian 
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